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1 .0  Introduction 


1.1  Purpose 

In  this  work  a  structural  theory  is  devised  for  the  analysis  of  laminated  anisotropic 
damped  plates.  The  plate  consists  of  laminated  face  sheets  sandwiching  a  single  damping  layer. 
The  displacement  degrees  of  freedom  used  to  construct  the  analytical  model  are  a  generaliza¬ 
tion  of  those  found  in  constrained  layer  theory.  Variational  principles  are  used  to  form  the 
equations  of  motion  which  are  solved  for  the  steady  state  excitation  of  simply  supported  semi¬ 
infinite  plates.  The  purpose  for  developing  such  an  analytical  model  is  to  provide  a  tool  for 
examining  the  effects  of  stress  coupling  and  compliant  layering  in  improving  the  performance 
of  damped  structures  loaded  in  flexure. 

1.2  Background 

Highly  damped  structure  is  built  by  combining  stiffness  and  damping  materials  into  a  sin¬ 
gle  structural  form.  An  example  of  such  a  composite  construction  is  a  planar  lamination  of 
stiffness  and  damping  layers.  Under  bending,  the  deformation  of  the  stiffness  layers  shear  the 
damping  layers  which,  because  of  their  viscosity,  dissipate  part  of  the  strain  energy.  Thus, 
resonant  vibrational  responses  are  controlled  by  the  conversion  of  vibrational  energies  into 
heat.  Damped  structures  based  on  this  principle  include  constrained  layer  treatments,  sandwich 
constructions  and  damped  plates  [1,2].  In  vibrational  environments  these  composite  components 
have  operational  characteristics,  durability,  and  service  life  that  far  exceed  that  of  components 
based  solely  on  stiffness  materials  [3,4]. 

Any  design  approach  that  increases  the  rate  or  amount  of  straining  in  the  damping  materi¬ 
als  of  a  composite  construction  has  the  potential  of  improving  the  structural  damping.  Towards 
this  end,  research  has  been  devoted  to  examining  the  use  of  discontinuous  stiffness  layers  in 
order  to  introduce  favorable  edge  effects  [5].  Also,  designs  have  been  investigated  that  increase 
the  thickness  deformation  of  the  damping  layers  [6,7].  Additional  design  approaches  can  be 
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found  in  the  use  of  anisotropy  and  lamination. 

When  anisotropic  stiffness  materials  arc  incorporated  into  damped  composite  construction, 
a  coupling  between  normal  and  shear  effects  (stress  coupling)  is  created.  Stress  coupling  pro¬ 
duces  components  of  strain  in  the  damping  layers  that  do  not  occur  in  isotropic  designs  while 
retaining  and  adding  to  those  components  that  already  do  occur.  Conceivably  this  increased 
straining  of  the  damping  material  would  lead  to  higher  energy  dissipation  and  improved  struc¬ 
tural  performance. 

Lamination  can  be  used  to  favorably  change  the  stiffnesses  of  a  structure.  For  instance, 
consider  a  damped  composite  plate  system  with  layered  stiffness  face  sheets  and  a  damping 
core.  If  the  stiffness  layers  adjacent  to  the  core  are  replaced  by  a  material  more  compliant  than 
the  outer  stiffness  layers  then  the  construction  would  have  virtually  unchanged  bending 
stiffness  but  directly  reduced  in-plane  facial  extensional  stiffness.  Under  transverse  vibration, 
the  axially  compliant  face  sheets  can  increase  the  amount  of  relative  in-planc  motion  between 
the  face  sheets,  thereby  leading  to  an  increase  in  the  rate  of  core  straining  and  energy  dissipa¬ 
tion. 

The  pioneering  work  in  the  use  of  stress  coupling  to  assist  the  control  of  dynamic 
responses  examined  the  performance  of  sandwich  constructions  [8],  The  results  of  this  study 
were  limited  since  the  analytical  model  was  based  on  displacement  fields  that  were  too  elemen¬ 
tary  to  fully  characterize  the  deformational  response  of  composite  constructions  with  compliant 
damping  layers.  This  shortcoming  was  rectified  in  a  later  work  which  examined  stress  coupling 
and  lamination  in  damped  bending  elements  [9],  An  additional  effort  in  this  area  studied  the 
use  of  stress  coupling  to  improve  the  damping  of  tension-compression  members  [10]. 

In  the  present  work  a  lamination  theory  is  formulated  that  is  applicable  to  a  general  class 
of  damped  plate  structures.  The  theory  contains  all  of  the  elements  necessary  for  examining  the 
effects  of  anisotropy  and  lamination  in  constrained  layer  treatments,  damped  sandwiches, 
damped  plates  and  other  damping  designs.  The  theory  also  includes  all  of  the  components  of 
translational  and  rotational  mass  inertia  plus  Timoshenko  shear  theory. 
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2.0  The  Theory  of  Anisotropic  Laminated  Damped  Plates 

2.1  Analytical  Model 

The  damped  plate  structure  consists  of  top  and  bottom  face  sheets  sandwiching  a  single 
damping  layer  (see  Figure  1).  The  face  sheets  are  layered  with  a  total  of  NT  layers  in  the  top 
face  sheet  and  NB  layers  in  the  bottom  face  sheet.  The  thicknesses  of  the  individual  layers  are 
designated  by  rj  for  the  top  layers,  tB  for  the  bottom  layers  and  tD  for  the  damping  layer. 
(Here  the  subscript  n  identifies  individual  stiffness  layers  while  the  superscripts  T  (top),  D 
(damping),  and  B  (bottom)  refer  to  specific  parts  of  the  structure).  The  right  handed  global 
coordinate  system  shown  in  Figure  1  and  used  in  the  development  consists  of  the  coordinates 
X\  and  x 2,  which  are  located  in  the  planar  mid-surface  of  the  damping  layer  (the  reference  sur¬ 
face),  and  the  transverse  coordinate  x3. 

To  analytically  model  the  described  structure  the  following  assumptions  arc  made: 

1.  The  in-plane  deformations  of  the  face  sheets  vary  linearly  through  the  face  sheet  thick¬ 
ness; 

2.  The  in-planc  deformations  of  the  damping  layer  vary  linearly  through  it’s  thickness; 

3.  The  in-plane  displacement  fields  are  continuous  across  the  interfaces  (perfect  bonding); 

4.  The  transverse  displacement  is  the  same  for  all  parts  of  the  cross  section. 

5.  The  moduli  of  all  of  the  materials  of  construction  can  be  treated  by  the  Complex 
Modulus  model; 

6.  The  material  model  for  the  stiffness  layers  is  transversely  isotropic  but  neglects  the  thick¬ 
ness  normal  stresses.  The  axis  of  isotropy  is  parallel  to  the  mid-surface; 

7.  The  material  model  for  the  damping  layer  is  isotropic  but  includes  only  shear  stresses. 


3 


Top  Layered  Face  Sheet 


NADC-90066-60 


l 


4 


Figure  1.  Damped  Plate  Structure 
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2.2  Theoretical  Development 

The  equations  of  motion  for  the  damped  plate  structure  are  derived  by  using  Hamilton’s 
Principle  in  conjunction  with  Reissner’s  Variational  Theorem.  Since  Hamilton’s  Principle  is 
only  applicable  to  conservative  systems,  the  material  properties  are  initially  treated  as  being 
purely  elastic  without  any  damping.  The  energy  integrals  are  then  formulated  and  minimized 
for  this  provisional,  fully  elastic  system.  Once  the  equations  of  motion  and  the  force- 
displacement  relations  are  derived,  damping  can  be  introduced  by  invoking  the  Correspondence 
Principle  in  which  the  elastic  moduli  are  replaced  by  the  complex  viscoelastic  moduli  of  the 
Complex  Modulus  model.  Application  of  the  damped  plate  model  is  therefore  limited  to 
steady  slate  harmonic  vibrations. 

Through  Hamilton’s  Principle  the  development  proceeds  by  finding  the  motion  that 
minimizes  the  integral 

*2 

l  =  j(K~R)dx  (2.2.1) 

where  K  is  the  kinetic  energy  of  the  system  and  R  is  Reissner’s  Functional.  The  integration  is 
over  the  time  variable  x  between  the  instants  of  time  x,  and  x2.  The  kinetic  energy  of  the  sys¬ 
tem  is  calculated  from 

*  =  yj \pUiiiidV  (2.2.2) 

where  p  is  the  mass  density,  V  is  the  volume  and  the  u,  arc  the  displacements.  (Throughout 
this  development  an  over  dot  will  indicate  differentiation  with  respect  to  time).  Reissner’s 
Functional  for  the  structure  is 

R  =  |(oy e,j -W (c,j ))d\;  -jT,u,dS  (2.2.3) 

where  and  c,j  arc  the  stress  and  strain  tensors,  W(0|;)  is  the  strain  energy  density 
expressed  in  terms  of  the  stresses,  the  T,  arc  the  tractions,  and  S,  is  that  portion  of  the  surface 
over  which  the  tractions  act. 
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The  integrands  of  equations  (2.2.2)  and  (2.2.3)  are  developed  from  the  assumptions  of  the 
analytical  model.  Starting  with  assumptions  1  through  4,  the  motion  of  the  structure  is 
expressed  in  terms  of  nine  unknown  displacement  degrees  of  freedom.  Using  these  functions 


the  displacement  fields  are  found  to  be: 

Top  Face  Sheet 

u1  =  u10(x1,x2.'t)-,-yfDaP(x1,x2,t)  +  (x3- YrD)a1r(x1,x2,'c)  (2.2.4) 

“2  =  “2  (X1,X2,T)+yt£)af(Xi,X2,T)  +  (X3-yrD)(xJ(x1,X2,T)  (2.2.5) 

u3  =  u°  (xx,x2,x)  (2.2.6) 

Damping  Layer 

u ,  =  u ,°  (x  i ,  x2,  T)  +  x3  aftx , ,  x2,  x)  (2.2.7) 

u2  =  u^(xl,x2,x)  +  x3aP(xi,x2,x)  (2.2.8) 

u3  =  u3(xl,x2,x)  (2.2.9) 

Bottom  Face  Sheet 

u\=ui  (*,-*2- x)~  y^Da1D(j:i,x2,T)+(x3  +  yrD)af(xI,A:2,x)  (2.2.10) 

U2  =  U2(xi,x2,'i)-^-tD  a?(x],x2,z)  +  (x3  +  jtD)a2(x],x2,x)  (2.2.11) 

«3  =  “3  (xux2,x)  (2.2.12) 


where  the  u,°  are  the  reference  surface  displacements,  the  aP  are  the  rotations  of  the  damping 
layer  about  the  reference  surface,  and  the  aj  and  a?  are  die  rotations  of  the  top  and  bottom 
face  sheets  (see  Figure  2). 

From  the  displacement  fields,  the  strains  in  the  reference  surface  and  the  plate  curvatures 
arc  computed  as 


e°  =  —  (u  0  +u  0  ) 
c«/  2  J  •* ' 
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Section  Parallel  to  the  x,  -  Coordinate  Line 


Section  Parallel  to  the  x2  -  Coordinate  Line 


Figure  2.  The  Displacement  Degrees  Of  Freedom. 
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*J=y(a£+a/,)  i,j= 1.2  (2.2.14) 

i,j =1,2  (2.2.15) 

^J-y(<-  +  aA)  ‘.7=1,2  (2.2.16) 

These  expressions  are  combined  to  yield  the  following  useful  quantities 

=  eS  +  yf°  (*£-*£)  i,j— 1,2  (2.2.17) 

*®  =  e<°-}'D  (*.?-**)  i,j=l,2  (2.2.18) 


The  strain  fields  are  now  computed  as 
Top  Face  Sheet 

tij=i[j+x3Kjj  i,j= 1,2  (2.2.19) 

e,3=“(a,T+«3°,)  <  =  1.2  (2.2.20) 

f-33  =  0  (2.2.21) 

Damping  Layer 

eij=etj+x3K-j  i,j= 1,2  (2.2.22) 

e13  =  y(a1D+«30,)  <  =  1.2  (2.2.23) 

e33  =  0  (2.2.24) 

Bottom  Face  Sheet 

Cij=*!j+x3K!}  ij =1,2  (2.2.25) 

e;3=y«x,fl  +  «30,)  <=1,2  (2.2.26) 

e33  =  0  (2.2.27) 


8 


NADC-90066-60 


Turning  to  assumption  6,  the  transversely  isotropic  constitutive  relations  for  the  stiffness 
layers  are  formed  in  the  local  layer  coordinate  system.  Tensorial  transformations  are  then  used 
to  reexpress  these  relations  in  the  global  coordinate  system.  The  matrix  form  of 

these  expressions  for  the  n  ’th  stiffness  layer  is 


<*11 

Clm 

^1122 

0 

0 

ClU2 

ell 

<*22 

C 1122 

C2222 

0 

0 

C2212 

£22 

o2"3 

= 

0 

0 

C  2323 

^2312 

0 

2^23 

<J,"3 

0 

0 

fn 

2312 

^1212 

0 

2Ei3 

of2 

C1112 

C22!2 

0 

0 

c  1212 

2e12 

(2.2.28) 


where  the  o'}  are  the  layer  stresses  and  the  C*jld  are  the  elastic  constants.  From  assumption  7, 
the  isotropic  constitutive  relation  for  the  damping  layer  is 


°23 

G  0  0‘ 

2^23 

°13 

= 

0  G  0 

2^13 

°12 

.0  0  G 

2Cl2 

(2.2.29) 


where  G  is  the  shear  modulus  of  the  damping  material. 

To  determine  the  layer  stress  resultants  the  constitutive  relations  are  integrated  over  the 
layer  thicknesses.  The  resulting  expressions  arc  then  inverted  to  yield  the  following  expres¬ 
sions 


Top  Layers 


Oi,‘  = - J— - —  + - - — ^ - —  x3 


'•j 


H 


nT 


H 


nT 


=  2 


nT 


nnT_fj3 
®i  3  t 


1  =  1.2 


Damping  Layer 


0,d2  = 


pD 

r  12 

.D 


12  w  D  _ 

dTT^i**3 


(tu) 


ij  =  1.2 


p D 
0,3  (D 


i  =  l,2 


(2.2.30) 

(2.2.31) 


(2.2.32) 


(2.2.33) 
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Bottom  Layers 


U nB  rr  nB  it  nB  1/ nB 
nB_  H  3  Fij  "H2  Mtj 

"  "  H- 


+ 


-H  fF-f  +  H  fMjf 


HnB 


x3 


°«3  - 


z?nB 

*i3 


1=1,2 


i,y  =  1,2 


(2.2.34) 

(2.2.35) 


where  the  F,"7,  F®,  and  F*5  are  the  layer  force  stress  resultants,  the  A/,"7,  Aff2,  and  A/*5  are 
the  layer  moment  stress  resultants  and  the  H"7 and  7//*®  are  geometric  quantities  computed  from 


Hf 

(2.2.36) 

Hf 

=  ytf  f(C.  +rf) 

(2.2.37) 

H? 

=  jWf((/If+1)2+/inr/iJ+1 

+(/«J)2) 

(2.2.38) 

HnT 

—  H  jj  nT  n  nT  it  nT 
— 17  j  It  3  “ 

(2.2.39) 

Hf 

(2.2.40) 

Hf 

=  jHf(h*+l+h*) 

(2.2.41) 

Hf 

=  l//f((/lf+])2  +  /if/if+1 

+(/jf)2) 

(2.2.42) 

HhB 

_ jj  nB  it  nB  jj  nB  jj  nB 

-n  |  n  3  —n  2  n2 

(2.2.43) 

in  which  the  hj  and  hj+l  are  the  x3  coordinates  of  the  bottom  and  top  surfaces  of  the  n  ’th  top 
layer  and  the  /if  and  /if+1  are  the  x3  coordinates  of  the  top  and  bottom  surfaces  of  the  n  'th 
bottom  layer. 

The  integrands  of  equations  (2.2.2)  and  (2.2.3)  are  now  computed  using  the  expressions 
developed  for  the  field  variables.  The  stress  resultants  are  included  by  performing  the  thickness 
integration  of  these  integrals.  The  resulting  expressions  arc  then  substituted  into  equation 
(2.2.1).  Taking  the  variation  of  I  with  respect  to  the  generalized  displacements  and  forces  and 
setting  the  coefficients  of  like  variations  to  zero  yields  the  governing  system  of  differential 
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equations.  These  equations  include  the  equations  of  motion,  the  force-displacement  relations, 
and  the  boundary  conditions  (the  initial  conditions  are  irrelevant  for  the  steady  state  problem). 
Following  this  procedure  the  force-displacement  equations  are  found  to  become 


r«T _ s^nT  / ij nT -mT  .  vj nT is-T \ 

Fij  - ^ijkl  \h  1  ^ki^H2Kkl) 

i,j=\,2,k,l  =  1,2 

(2.2.44) 

isnT—.r'nT  ( Li nT -pT  i  nT msT \ 

Mij  -LijkiyH  2  Zkl+H  3 

i,j=\,2,k,l  =  1,2 

(2.2.45) 

F”f,  =  C$(Hfz%+H?K%) 

i,j=\,2;k,l  =  1,2 

(2.2.46) 

i,j=\,2,k,l  =  1,2 

(2.2.47) 

Fr? =t?Ciy3  («/ +u3j) 

(,>=1,2 

(2.2.48) 

**  ='?C!fj3  (a /+«&) 

(,>=1,2 

(2.2.49) 

F?2  =2G  ^e,0. 

(2.2.50) 

^f2  =  ^G  (r°  )3/ff2 

(2.2.51) 

F^=Grw(a1D  +  «3°1) 

/  =  1,2 

(2.2.52) 

Summing  these  expressions  over  the  top  and  bottom  face  sheets  leads  to 
Top  Face  Sheet 

Ffj=Aijuukj  +  ^tD AljuCtij  +(-yfD AjjU  +BjjU)al j  i,j= 1,2  ;  k,l  =  1,2  (2.2.53) 

Fl ,  =  EjjiaJ  +  u?j)  i,j= 1,2  (2.2.54) 

+  ^tDBljkla?j  +  (-±tDB?jkJ  +D[jkl)al,  i,j= 1,2  ;  k,l  =  1,2  (2.2.55) 

Bottom  Face  Sheet 

+  *'•>=!• 2  ;  k’l  =  ]’2  (2.2.56) 

F»i  =  Ejj(  a/  +  u3°i/)  i,  >=1,2  (2.2.57) 

M*  =  -  j ®*/«U  +  BfjU  +  D,5u)a£,  i , >=1 , 2  ;  A ,  /  =  1 , 2  (2.2.58) 
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where  the  F?,  Ffj,  Mj}  and  M?  are  the  face  sheet  force  and  moment  stress  resultants  and 


Ail- ECjl  H? 

n  =  1 

i,j=\,2\  k,l  =  1,2 

(2.2.59) 

N 

d  V  _  r,  s*nT  jjnT 

aijU  ~  L  ijkl  H  2 
*  =  1 

1,7=1, 2;  A,/  =  1,2 

(2.2.60) 

N 

n T  _  V1  /'"I" 

h(-ijklt1  3 
*  =  1 

1, 7=1, 2;  k, 1  =  1,2 

(2.2.61) 

hj 

II 

r> 

i',7=l,  2 

(2.2.62) 

M 

A  B  _  V*  /""fl  r/nB 

A1JU  ~  l*  C‘jU  n  1 

n  —  1 

i,j=l,2;  k, 1  =  1,2 

(2.2.63) 

M 

D B  V*  iin# 

Lijkl  n  2 

«  =  1 

1,7=1, 2;  k, 1  =  1,2 

(2.2.64) 

M 

n#  __  v1  rjmB 

L>ijki~  1a  c  ijkl  H  3 

n  =  1 

i,7=l, 2;  A: , /  =  1, 2 

(2.2.65) 

Efj=  ZZchjs 

n  =  1 

1,7=1, 2 

(2.2.66) 

Using  equations  2.2.53  to  2.2.58  in  the  variational  relations  leads  to  the  following  form 
of  the  equations  of  motion 


-FU\;-fL-^2.2-^1+M«i0+/?aiD+/[ar+/iaf  =  0 

7=1,2 

(2.2.67) 

~F2jj~F2jJ-F?Z,2  -P 2+M“2  +/f  of  +  /[d2r+/f  df  =  0 

7=1,2 

(2.2.68) 

-^3,  -  ^3,  -  F?x<  -  Py +  ^«3°  =  0  1  =  1 ,2 

(2.2.69) 

(-F \j j  +  F\j j)  +  Ftf  —Myi,2  F l^ii\  +  /2>  ctp  +  f 2  o^+/f  df  =  0  7=1,2  (2.2.70) 
\‘D(-FljJ+FB2Jj)  +  F&-M?2']  +/f«2°+/?df  +  /[dr+/fd2fl  =  0  7=1,2  (2.2.71) 
~tD  F^j+Fl 3  -Mfjj  +  lW  +lld?+ll*i=0  7=1,2  (2.2.72) 
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F2jj+F23  ~M2jj  +  1]>*2  +I2&2 +1I&2 =0  7=1*2  (2.2.73) 

-\*D  ^j+Ff3-<;j+/f«i°+/faiD+/faf=0  7=1,2  (2.2.74) 

-y'D  f2jj+f23  -Wf;j  +/?«20+/f  af+/f  ^2fl=0  7=1,2  (2.2.75) 

where  the  P,  are  the  applied  tractions  and  the  inertial  constants  are  computed  from 


mt=i,p  Itl 

n  =  1 

(2.2.76) 

A/B=Xp»'-f 
«  =  1 

(2.2.77) 

M  =Mt  +pDtD  +Mb 

(2.2.78) 

,T  y  T  T(  hn+\  +hn  1  Ds 

•  1  -  Zj  tn  Pn  (  -  ~l  ) 

n  =  1  z  z 

(2.2.79) 

/f  =lrD(Mr-MB) 

(2.2.80) 

M  hB  ,  +/lB  1 

/?=l Mr"* n  +\tD) 

«=i  z  z 

(2.2.81) 

7r  _  1  fD,r 
'2  - y'  *  1 

(2.2.82) 

/£  =  (*  tD)2(MT  +MB)  +  pD  (f°)3 

2  2  P  12 

(2.2.83) 

/5~f0/? 

(2.2.84) 

'3=1  <(-|'D»Ip.V-'DP.IHf +P.rHf) 
n  =  1  Z 

(2.2.85) 

'!  =  S  ((V/p^  +  'V^  +  P>f  ) 

M  =  1  ^ 

(2.2.86) 

where  pj  and  pB  are  the  layer  mass  densities. 
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The  natural  boundary  conditions  are  found  to  be 


(F\jnj  +F?2n2+FBljnj)6u ,°  =0  >= 1,2 

(2.2.87) 

(Fljfij  +F?2nx+Fljnj)8u%  =0  >= 1,2 

(2.2.88) 

(Fiini  +Fnnl  +  F%3n2+ F?3nt )  5 u  3°  i  =  1 ,2 

(2.2.89) 

(|rDF[>n>+Mf2«2-|tDFf/n7)6aP=0  >= 1,2 

(2.2.90) 

(jrDF^n7+Mf2rti-|rDFf7ny)5a2D=0  >= 1,2 

(2.2.91) 

(~j‘D  F\jrij  +A/[,«;)8a,r=0  >= 1,2 

(2.2.92) 

+MT2jnj)boJ =0  >=1,2 

(2.2.93) 

(|rDFfyny+M?7ny)5af=0  >=1,2 

(2.2.94) 

(i-fDF«n7+Aff7n;)5af=0  >=1,2 

(2.2.95) 

where  6  is  the  variational  symbol  and  the  rij  are  the  direction  cosines  of  the  normals  on  the 
structural  boundaries. 

At  this  point  the  force-displacement  relations  are  substituted  into  the  equations  of  motion. 
This  yields  a  set  of  nine  displacement-equilibrium  equations  the  unknowns  of  which  are  the 
nine  functional  displacement  degrees  of  freedom.  Solutions  to  specific  problems  are  found  by 
applying  the  appropriate  set  of  boundary  conditions  and  solving  these  displacement-equilibrium 
equations. 
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3.0  Semi-Infinite  Simply  Supported  Plates 


3.1  Structural  Description 

The  equations  of  the  anisotropic  laminated  damped  plate  theory  are  extremely  complex 
and  are  not  readily  solved  by  closed  form  or  series  solution  methods.  For  instance,  individual 
equations  of  the  theory  contain  odd  and  even  derivatives  of  the  same  spatial  variable,  which 
precludes  a  Fourier  series  solution.  Semi-infinite  simply  supported  plates,  however,  are  an 
important  class  of  problems  for  which  a  Fourier  series  solution  can  be  found  and  in  which 
stress  coupling  and  lamination  can  be  studied. 

Consider  z  plate  of  finite  length  a  in  the  Xj  direction  and  of  infinite  extent  in  the  x2 
direction.  On  the  xt  =  0  and  x^=a  edges  the  plate  is  taken  to  be  simply  supported.  Further¬ 
more  the  applied  loads  on  the  plate  are  limited  to  be  uniform  in  the  x2  coordinate.  For  such  a 
system  all  of  the  field  variables  are  independent  of  the  x2  coordinate,  so  that  all  of  the 
differentiations  with  respect  to  this  coordinate  vanish.  This  simplifies  the  equations  of  motion 
which  in  matrix  notation  reduce  to 

[M][U)  +  [D)[u]  =  [P]  (3.1.1) 

where  [M]  is  the  mass  matrix,  [D]  is  a  differential  operator  matrix,  [u]  is  a  vector  of  unknown 
displacement  functions  and  [ P ]  is  a  load  vector.  The  elements  of  these  matrices  arc  given  in 
Appendix  A. 

The  natural  boundary  conditions  for  this  structure  reduce  to  the  following  edge  conditions 
at  *1=0  and  x,=a 


Ffi=0 

(3.1.2) 

F?,=0 

(3.1.3) 

FTm  +F?2  +Ff2=0 

(3.1.4) 

u3°=  0 

(3.1.5) 
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|rDF[2+Mf2-|rDFf2=0  (3.1.6) 

M\l=0  (3.1.7) 

MBn  =0  (3.1.8) 

-\tDpT^  +Mn  =0  (3.1.9) 

|/DF?2+M?2=0  (3.1.10) 

3.2  Solution 

Equation  (3.1.1)  is  solvable  by  the  Fourier  series  method.  To  apply  this  method  assume 
the  following  series  expansions  for  the  displacement  degrees  of  freedom 

u,°=  £f/?cos(^^)e,iiT  (3.2.1) 

m= 1  a 

U2°  =  £  l/?  cos(^l)*-  (3.2.2) 

m=\  & 

«30  =  £  U"  sm(^±)eiQz  (3.2.3) 

m-l 

aiD=  £AT°  cos(^±)e‘‘*  (3.2.4) 

m=l 

aD=  £^mo  cos(-^^-)e‘at  (3.2.5) 

m=l  ^ 

atT=  £^i,rcos(^^-)e‘nT  (3.2.6) 

m=l  ^ 

0.2=  £/l^rcos(-^i)e‘at  (3.2.7) 

m=l  ^ 

af=  £4^  cos(-^^-)e'nT  (3.2.8) 

m=l  ^ 
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af=£A?  cos(-^V<* 

m- 1  ^ 


(3.2.9) 


where  the  superscripted  constants  are  Fourier  coefficients  and  £2  is  the  frequency  of  the  steady 
state  excitation.  It  is  easily  shown  that  these  expressions  satisfy  the  boundary  conditions. 

The  harmonically  varying  excitations  (with  respect  to  time)  are  also  expressed  in  terms  of 
Fourier  series  expansions 


m  =  l 


m- 1 


m- 1 


mKX] 


-)eiDx 

(3.2.10) 

(3.2.11) 

)eiC* 

(3.2.12) 

where  the  P,m  are  the  Fourier  coefficients  determined  from  the  Fourier  formulacs. 

Substituting  the  above  expansions  into  equation  (3.1.1)  results  in  an  infinite  number  of 
uncoupled  equations  that  can  be  grouped  into  sets  by  common  indicial  values.  Thus  a  set  of 
nine  equations  and  nine  unknowns  is  obtained  for  each  indicial  value  where  the  unknowns  of 
these  equations  are  the  Fourier  coefficients  of  the  displacement  scries.  Expressing  these  equa¬ 
tions  in  matrix  form  leads  to  the  following  general  expression  for  each  indicial  value 

-tf[M  ]  [Um  ]  +  [Bm  }{Um}  =  \Pm  ]  (3.2.13) 


where  [Um]  is  a  vector  of  Fourier  displacement  coefficients,  [ Bm  ]  is  a  modal  stiffness  matrix 
whose  elements  are  determined  by  the  material  and  geometric  properties  of  the  structure,  and 
[Pm]  is  a  vector  of  the  Fourier  loading  coefficients.  The  elements  of  these  matrices  are  given 
in  Appendix  A. 

The  viscoelasticity  of  the  original  structure  may  now  be  introduced  by  invoking  the 
Correspondence  Principle.  This  is  done  by  replacing  the  clastic  material  moduli  with  complex 
viscoelastic  moduli.  The  clastic  field  variables  arc  then  reinterpreted  as  complex  harmonic  vari¬ 
ables  and,  as  assumed  in  the  Fourier  scries  expansion  of  the  loadings,  the  excitation  is  limited 
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to  steady  state  harmonic. 

3.3  Analytical  Approach 

The  analysis  can  be  completed  in  several  ways  depending  upon  the  type  of  information 
desired.  For  instance,  the  dynamic  response  of  a  damped  plate  to  a  specific  excitation  can  be 
found  through  the  direct  solution  of  equations  (3.2.13).  The  analysis  begins  by  decomposing 
the  spatial  portion  of  the  excitation  into  it's  Fourier  components.  If  necessary  a  suitable  level 
of  truncation  is  determined.  Equations  (3.2.13)  are  then  formed  and  solved  for  each  applicable 
indice.  The  displacement  degrees  of  freedom  are  then  synthesized  from  the  Fourier  expansions. 
Once  the  displacements  are  found,  then  any  field  variable  of  interest  can  be  computed.  If 
response  spectrums  are  desired,  the  steps  of  this  method  are  repeated  throughout  the  frequency 
range  of  interest.  Frequency  dependent  material  properties  arc  accounted  for  by  updating  their 
values  at  each  step  in  the  analysis. 

If  however,  the  modal  loss  factors  arc  to  be  determined  then  the  Damped  Forced  Vibra¬ 
tion  Mode  Method  [11]  is  applied.  In  general,  high  system  damping  is  nonproportional  so  that 
there  is  a  coupling  between  undamped  normal  modes,  but  a  unique  value  of  modal  damping 
can  still  be  found  by  considering  a  special  case  of  forced  vibration.  In  particular,  consider  har¬ 
monic  loads  that  are  proportional  to  the  local  mass  inertia  forces  and  that  are  in  phase  with  the 
local  velocity.  Then 

[Pm]  =  -iT\Cl2[M][Um]  (3.3.1) 

For  this  loading  Equation  (3.2.13)  becomes 

[flm][Cn-«2(l+n)[M][f/m]  =  [0]  (3.3.2) 

which  can  be  rearranged  into  the  form  of  a  complex  eigenvalue  problem 

[[flm]-Q2(l+/Ti)fM]][(/'"]  =  [0]  (3.3.3) 

From  this  last  equation  the  natural  frequencies,  loss  factors,  and  complex  mode  shapes  arc 


extracted. 
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4.0  Applications 


4.1  Structural  Description 

The  semi-infinite  damped  plate  examined  in  this  analytical  study  has  a  width  of  25.4  cm. 
The  top  and  bottom  face  sheets  of  the  plate  consist  of  6  stiffness  layers  with  each  layer  having 
a  thickness  of  0.1/25  mm.  The  damping  layer  has  a  thickness  of  .0965  mm.  The  stiffness 
layers  consist  of  IM6-3501  carbon-epoxy  with  a  fiber  volume  fraction  of  60%.  The  properties 
of  this  material  are  shown  in  Table  1  where  the  frequency  dependence  of  the  properties  are 
neglected  since  the  overall  structural  damping  is  dominated  by  the  damping  layer.  The  damping 
layer  consists  of  Type  112  Scotchdamp  SJ2015x.  The  frequency  dependence  of  the  moduli  of 
this  material  arc  accounted  for  in  the  analysis.  A  plot  of  these  properties  is  shown  in  Figure  3. 
The  mass  density  of  the  damping  material  is  .98  gm/cc. 

To  study  the  effects  of  stress  coupling  and  compliant  layering  on  structural  damping,  the 
fiber  reinforced  layers  adjacent  to  the  damping  layer  arc  given  off-axis  orientations  with  respect 
to  the  Xj  coordinate  direction.  This  arrangement  of  layers  retains  static  stiffness  by  keeping  the 
fibers  in  the  outer  plys  oriented  along  the  plate  width.  However,  with  off-axis  orientations  for 
the  inner  layers,  the  top  and  bottom  face  sheets  can  individually  be  unbalanced  and  unsym- 
metrical  leading  to  full  coupling  effects  (extcnsional-shcar,  extcnsional-bcnding/twisting  and 
bending-twisting).  Also,  the  off-axis  orientation  makes  the  inner  layers  more  compliant  with 
respect  to  the  x,  coordinate  direction.  When  the  face  sheets  are  bonded  together  by  the  damp¬ 
ing  layer  the  overall  structural  coupling  is  enhanced  or  diminished  depending  upon  the  layer 
orientations  and  the  stiffness  of  the  damping  layer.  For  instance,  with  stiff  damping  layers  the 
Love-Kirchoff  hypothesis  is  valid  so  that  the  structural  coupling  of  the  damped  plate  is  identi¬ 
cal  to  that  predicted  by  lamination  theory.  However,  for  flexible  damping  materials  the  top  and 
bottom  face  sheets  arc  not  constrained  to  act  in  unison  so  that  more  complicated  coupling 
phenomena  can  occur. 
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Table  1.  Material  Properties  Of  IM6-3501  Carbon-Epoxy. 


Axial  Extensional  Modulus 
Transverse  Extensional  Modulus 
Axial  Poisson's  Ratio 
Axial  Shear  Modulus 
Transverse  Shear  Modulus 
Axial  Loss  Factor 
Transverse  Loss  Factor 
Shear  Loss  Factor 
Mass  Density 


148, GPa  (21.5  Msl) 

8.96  GPa  (1.30  Msl) 

.35 

4.48  GPa  (0.65  Msl) 

2.07  GPa  (0.30  Msl) 

.00128 

.0110 

.0110 

1.52  gm/cc  (.000142  #sec?/INA) 
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Real  &  Imaginary  Shear  MODULUS  (MPa) 
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Figure  3.  Material  Properties  Of  Scotchdamp  SJ  2015  X.Type  112  @  20°C. 
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The  notation  used  to  specify  the  structural  arrangement  of  the  damped  plates  is  identical 
to  that  used  for  laminations  of  advanced  composites  except  for  the  addition  of  the  symbol  d 
which  will  indicate  the  presence  of  a  damping  layer.  For  instance,  the  baseline  structure  for 
this  study,  so  called  because  it  does  not  include  stress  coupling  or  compliant  layering  effects,  is 
designated  0 6ld  /06. 

4.2  Numerical  Results 

The  first  analysis  will  examine  structures  with  layups  of  O4/02/  d  /02/O4  and 
O4/02/d/-02/O4  where  the  layer  orientation  angle  0  is  varied  between  0  and  90  degrees.  The 
Forced  Mode  Method  is  used  to  compute  the  loss  factors  of  the  first  three  bending  modes  of 
these  structures.  Plots  of  the  results  arc  shown  in  Figures  4  and  5.  The  plots  show  that  for  this 
structure  the  stress  coupling  and  compliant  layering  have  no  effect  on  the  loss  factor  of  the 
fundamental  mode  but  that  these  conditions  lead  to  increases  in  damping  in  the  two  higher 
modes.  Also,  it  is  seen  that  even  though  stress  coupling  vanishes  at  0=90  degrees,  the  highest 
loss  factors  are  achieved  at  this  orientation.  Thus  the  use  of  compliant  inner  layering  appears  to 
have  more  of  an  effect  on  damping  than  stress  coupling.  Therefore,  the  rest  of  the  analyses 
presented  here  will  consider  laminated  designs  without  stress  coupling.  In  addition  the  analyses 
will  also  be  limited  to  the  study  of  beams  since  the  second  in-plane  dimension  is  no  longer 
needed  in  the  absence  of  stress  coupled  lay-ups.  (Note  that  in  other  structural  designs  stress 
coupling  may  prove  to  have  a  greater  role  in  damping.  Also  note  that  most  of  the  increase  in 
damping  is  reached  by  0=60  degrees.  This  fact  may  have  importance  in  the  fabrication  of 
damped  structure.) 

Figure  6  shows  the  loss  factors  of  four  different  damped  beams  for  the  first  five  bending 
modes  of  vibration.  Here  it  is  seen  that  there  is  little  or  no  gain  in  damping  for  the  fundamen¬ 
tal  mode  but  that  in  the  higher  modes  the  compliant  inner  layer  laminates  have  significantly 
greater  loss  factors.  (The  matching  of  natural  frequencies  of  the  plates  indicates  that  the  gain 
in  damping  is  not  due  to  changing  material  properties). 
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Figure  4.  Modal  Damping  for  Bending  Modes  1  To  3. 
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The  ineffectiveness  of  the  compliant  layering  in  the  fundamental  mode  is  attributed  to  the 
particular  mix  of  design  variables.  Later  in  the  text  (Table  2),  a  design  is  shown  that  improves 
the  damping  in  the  fundamental  mode  but  only  by  resorting  to  commercially  unavailable  ply 
thicknesses.  In  other  design  situations  the  variables  (dimensions,  material  properties,  etc.)  may 
be  such  that  improvements  in  performance  can  be  found  across  the  spectrum. 

The  goal  of  a  damping  design  is  to  reduce  resonant  stresses  and  displacements.  This  is 
achieved  by  increasing  the  structural  loss  factor  which  in  the  laminated  designs  is  accomplished 
by  sacrificing  static  stiffness  (i.e.  through  the  use  of  90  degree  layer  orientations).  It  is  neces¬ 
sary  then  to  verify  that  the  structural  response  actually  decreases  in  the  highly  damped  but 
more  flexible  laminated  designs.  To  analytically  test  the  response,  the  structures  arc  subjected 
to  forcing  functions  that  approximately  excite  the  resonant  response  (the  approximation  is  intro¬ 
duced  by  not  accounting  for  the  negligible  in-plane  and  moment  components  of  the  load  vector 
that  are  required  by  the  Forced  Mode  method  for  a  strict  proportionality  to  the  inenia  loading). 
Figure  7  shows  the  result  of  this  computation  where  the  amplitude  of  the  transverse  displace¬ 
ments  have  been  normalized  with  respect  to  the  modal  response  of  the  baseline  plate.  Except 
for  the  fundamental  mode  where  virtually  no  improvement  is  achieved,  the  analysis  predicts 
reduced  resonant  responses. 

The  optimal  design  for  vibration  resistance  is  determined  through  a  balance  of  inenial, 
stiffness  and  damping  properties.  For  a  given  modal  excitation,  the  baseline  structure  has  a 
unique  optimal  configuration.  But,  since  the  damping  material  properties  arc  frequency  depen¬ 
dent,  the  optimal  design  at  one  mode  of  response  docs  not  necessarily  imply  an  efficient  struc¬ 
ture  over  the  rest  of  the  spectrum.  In  making  comparisons  between  the  response  of  the  lam¬ 
inated  designs  and  of  the  baseline  structure,  there  is  the  concern  that  the  baseline  structure  is 
not  in  an  optimal  modal  configuration.  To  test  this  concern,  analysis  is  performed  to  establish 
the  optimal  baseline  design  for  each  mode  of  response.  (Note  that  the  ability  to  actually  build 
an  optimal  design  is  limited  by  the  commercial  availability  of  ply  thicknesses.  The  analysis  is 
performed  here  for  the  academic  reason  of  establishing  the  efficacy  of  the  new  design 
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Figure  6.  Modal  Damping  For  Bending  Modes  1  To  5. 


NADC-90066-60 


o  ^ 

i  - 

?  s 

cr 

o  l 
§  -S 

^  Q. 

E 

_  to 

g  e 

O  (0 


uoj^aaijaQ  paojoj  pazijBUJJO^ 


u. 


27 


igure  7.  Forced  Center  Deflection  For  Bending  Modes  1  To  5. 
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approach.)  The  overall  thickness  of  the  plate  remains  fixed  at  2.17  mm.  Given  this  constraint, 
the  thicknesses  of  the  damping  and  stiffness  layers  are  varied  to  determine  the  optimal 
configuration.  Once  an  optimal  design  is  determined  for  the  baseline  structure,  the  compliant 
layered  designs  are  optimized  by  keeping  the  thickness  of  the  damping  layer  constant  while 
varying  the  thicknesses  of  the  0  degree  layers  and  90  degree  layers.  The  results  of  this  analysis 
are  shown  in  Table  2  where  it  is  seen  that  for  each  mode  of  response  a  laminated  design  is 
found  that  is  superior  to  the  baseline  design  in  increasing  the  loss  factor  and  reducing  the 
resonant  response. 

The  controlling  parameter  in  increasing  the  damping  in  the  compliant  layered  designs  is 
the  extensional  modulus  of  the  compliant  layers.  This  is  seen  in  Figure  8  where  the  modulus  of 
the  inner  layers  is  varied  parametrically  as  a  percentage  of  the  modulus  of  the  outer  layers.  The 
loss  factor  directly  increases  with  decreasing  modulus.  This  modulus  also  controls  the  phase 
lag  between  the  damping  layer  rotation  (af)  and  the  other  displacement  degrees  of  freedom 
(which  respond  approximatly  in-phase).  Figure  9  shows  that  this  phase  lag  increases  with 
decreasing  modulus.  This  means  that  as  the  beam  passes  through  it’s  static  equilibrium  posi¬ 
tion,  the  face  sheets  arc  increasingly  displaced  longitudinally  with  respect  to  one  another  as  the 
compliant  layer  modulus  decreases. 
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Table  2.  Optimized  Modal  Designs 


Mode 

Type  (1) 

td/t 

/'// 

t2lt 

Loss  Factor 

Displacement  (2) 

1 

B 

6.5 

0.0 

93.5 

.311 

1.00 

C 

6.5 

18.5 

75.0 

.341 

.80 

2 

B 

3.2 

0.0 

96.8 

.296 

1.00 

C 

3.2 

23.1 

73.7 

.332 

.76 

3 

B 

2.3 

.294 

1.00 

C 

2.3 

27.7 

70.0 

.337 

.72 

4 

B 

1.4 

0.0 

98.6 

.263 

1.00 

C 

1.4 

23.1 

75.5 

.300 

.77 

(1)  Structural  design; 

B  =  Optimized  baseline  design 
C  =  Optimized  compliant  layer  design 

(2)  The  transverse  displacement  is  normalized  with  respect  to  the  modal  response  of  the  base 
line  plate 
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Figure  8.  Modal  Damping  Vs  The  Modulus  Of  The  Compliant  Layer. 
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5.0  Conclusions 

In  order  to  examine  the  use  of  anisotropy  and  lamination  in  damped  plates  a  structural 
theory  was  developed  and  applied  to  a  simple  but  representative  structural  system.  The  analyti¬ 
cal  study  revealed  that  compliant  layering  can  increase  the  efficiency  of  damping  designs. 

Compliant  layering,  which  is  the  replacement  of  face  sheet  material  with  a  less  stiff 
material  at  the  interface  of  the  face  sheets  and  the  damping  layer,  affects  the  dynamic  response 
of  the  plate  through  the  alteration  of  stiffness  properties.  Compliant  layering  reduces  the  face 
sheet  axial  stiffness  without  altering  the  overall  structural  bending  stiffness.  This  creates  a 
mechanism  for  increasing  the  rate  of  shearing  in  the  damping  material  by  increasing  the  rela¬ 
tive  longitudinal  displacements  of  the  face  sheets.  The  shearing  rate  and  the  associated  energy 
dissipation  were  found  to  increase  as  the  modulus  of  the  compliant  layer  was  reduced.  How¬ 
ever,  there  is  a  limitation  to  this  process  since  the  stiffness  of  the  compliant  layer  must  be  high 
enough  to  confine  the  shear  deformation  to  the  damping  layer. 

Comp)’ ant  layering  can  also  be  used  to  reduce  the  weight  of  damped  structures  since 
compliant  materials  are  generally  less  massive  than  stiff  materials.  For  instance,  metallic  face 
sheets  that  incorporate  a  glass-epoxy  compliant  layer  can  have  improved  dynamic  resistance  at 
a  reduction  in  weight.  This  same  effect  can  be  achieved  by  merely  removing  some  of  the 
material  on  the  inner  side  of  the  face  sheets  through  grooving,  waffling  or  scoring  this  surface. 

Compliant  layering  introduces  challenges  to  the  fabrication  process  since  it  involves  cither 
the  mating  of  dissimilar  materials  or  the  unbalancing  of  quasi-isotropic  laminates.  Also,  there 
will  be  additional  steps  in  the  laminate  fabrication  which  will  add  to  the  cost  of  building  these 
components.  Nevertheless,  depending  upon  the  total  cost  of  construction,  compliant  layering 
offers  an  important  design  option  in  the  use  of  damped  bending  elements. 
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Appendix  A  -  Matrix  Elements 


The  elements  of  the  mass  matrix  [M]  arc 


[M]= 


M  0 
0  M 
0  0 
/?  0 
0  /? 
/[  o 
0  /[ 
/?  0 
0  /? 


0  /?  0  /[  0  /? 

0  0  /?  0  /[  0 

M  0  0  0  0  0 

0  /?  0  l\  0  /? 

0  0  /?  0  0 

0  /£  0  /[  0  0 

0  0  /£  0  /£  0 

0  /f  0  0  0  /f 

0  0  /f  0  0  0 


0 

/? 

0 

0 

/* 

'2 

0 

0 

0 

/« 

'  3 


The  elements  of  the  operator  matrix  [D  ]  are 


Olld.l, 
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where  the  DtJ  are 


^  11 = — (ATih  +^fni) 

D\2  =  D2\  = -(A  1112  +  A?n2) 

D  \%=D  ix  =  —tD  (~ATmi  +^fii2) 

D tf=D2x  =  (—tD A]u2  -Bjui) 

Di9=D9i=-(—tDAfn2  +£?m ) 

Z)  22  =  “04  1212  +i4f2!2  +GtD) 

D2A~B>A2~  (“A  i"ii2  +^U12) 

£,26  =  ^62=('J,D^1112  “^[112) 

Dn~D%2=-(,—iD  A^X\2  +  £?ii2) 

D33  =  -(£n  +£fi  +GtD) 

D-S4  =  -DAJ  =  -G  tD 

£,37  =  -£>73  =  -£i2 

Dy)  =  -D92  =  -E\2 

Dmx  =  -(—Id)  (A[U1  +i4?ul) 

DA5  =  Dsa  =  -(-^(D)  (A[h2 +Afu2) 

D47  =  D74=^-r0(^-r£>y4[u2  -£[112) 

£>49=1)94=  -^tD(-^tDAiiX2  +Bun) 

n  _  (—tD\(AT  +Ab  )  G(,D )3 

£>551  ~~'~2  '  1212  +  A  1212)  J2 

£>56  =  £>65  =  -^tD(^tDA  1112  -£[112) 

Ds%  =  Di5  =  —tD  (—tD  Axu2  +  B\\\2) 
D(A]  =  (-^tD)  AT\u\  +{D BJ]U  -D nil 


B>u=D4l  =  ~tD(-Ajni  +  A  fin) 
B>\6=D6X  =(—tD  A]m  -£nn) 

£,18  =  ^8l“-('^,C'4?lll  +£?lll) 

£>25  =  £>52=  y>D(_4[2 12  +A?212) 

D27  =  £>72=(‘J>°'4[212  “£[212) 

D2i)  =  D<)2  =  -(^tD  AB2\2  +£? 212  ) 

£>36  =  ~£>63  =  _£n 
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D  552  =  Gt^ 

Dst  =  Dis=  -^tD  (  —  tD  A]2\2  ~BT\ 212  ) 
£>59  =  £>95  =  {—1°  AB2\2  +B\ 212) 

£>  662  =  £  II 
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i  2 

E  611  =^761  =-(ytD)  Af112  +t°B fn2  “£>1112 

D  672  =  £>762  =  ^12 

Dn\=-(~2tD)  ^1212  +tD  B\ 212  -£>1212 

D  772  =  £  22 

Da\=-{-^tD)  i4fm  -r^fifin  -£>fm 

£)882  =  £?1 

£)*91=£)981=_(‘2,D)  ^ fl  12  -jfl^fll2  -^1112 

Dg92  =  Dgg2  =  Ei2 

£)991=-(y,C)  ^?212  -tC^?212  ~£>?212 

D992  =  E22 

and  where  dyl  and  d,n  denote  the  first  and  second  order  differentiation  with  respect  to  the  spa¬ 
tial  variable  xu 
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The  displacement  degree  of  freedom  vector  is 


The  load  vector  is 

>1' 

Pi 

Pi 

0 

[/>]=  0 
0 
0 
0 

.0. 

The  Fourier  displacement  coefficient  vector  is 
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The  Fourier  loading  coefficient  vector  is 
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The  elements  of  the  symmetric  [flm]  matrix  are 
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where  the  BX  are  computed  from  either 

ij*  3 

or 

-B^B^D^lf) 

with  the  exception  of  the  following  elements 
B^ 

a 

B?=-DMX(—j +D„2 
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b*=-d5SI(^)2+d552 
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